UdeM-GPP-TH-09-179 



CP Violation in Three-Body Chargino Decays 

G\ 

Makiko Nagashima/'Lj Ken Kiers, 2 '[] Alejandro Szynkman, 1 '^ 
^ David London, 1 '^ Jenna Hanchey, 2 '^ and Kevin Little 2 'p*1 

J> 1 Physique des Particules, Universite de Montreal, 

CP. 6128, succ. centre-ville, Montreal, QC, Canada H3C 3J7 
2 Physics and Engineering Department, Taylor University, 
236 West Reade Ave., Upland, IN 46989, USA 
(Dated: November 11, 2009) 

^ Abstract 

CP violation in supersymmetry can give rise to rate asymmetries in the decays of su- 
Qh persymmetric particles. In this work we compute the rate asymmetries for the three-body 

chargino decays xt ~> X\HH, xf -> xfZZ, xt -> xfW + W~ and xt -> XiZH. Each 
of the decays contains contributions mediated by neutral Higgs bosons that can possibly 
£Sj go on shell. Such contributions receive a resonant enhancement; furthermore, the strong 

phases required for the CP asymmetries come from the widths of the exchanged Higgs 
bosons. Our results indicate that the rate asymmetries can be relatively large in some 
cases, while still respecting a number of important low-energy bounds such as those com- 
ing from B meson observables and electric dipole moments. For the parameters that we 
consider, rate asymmetries of order 10% are possible in some cases. 
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Supersymmetry (SUSY) has been proposed as a solution to the hierarchy problem 
in the standard model (SM). SUSY is widely thought to be the physics that lies 
beyond the SM, and it is hoped that it will be discovered in the future at the LHC 
or at a linear collider. 

In SUSY theories each ordinary fermion and gauge boson has a superpartner, 
respectively of spin and spin |. SUSY also includes a charged Higgs boson, the 
H . The W~ and H~ each have a fermionic partner, known as charginos. These two 
charginos can mix, resulting in two mass eigenstates Xi an d X2 whose masses can 
be very different. Here we adopt the standard notation m^- > rn--. In this paper, 

we study CP violation in the decay xt -> xfXY, where XY = HH, ZZ, W+W~ 
or ZH (H is a Higgs boson). This analysis is an extension of the work performed 
in Refs [JJ [2J. Throughout this paper we assume that the CP-conserving SUSY 
parameters are known, but that those which violate CP remain to be measured. 
Some previous studies of CP violation involving charginos may be found in Refs. [21 

n ei Ei mm ei Emm]. 

CP violation in SUSY has been studied extensively at low energies, in meson 
mixing [12] . in the 5-meson system [13J and in electric dipole moments (EDMs) 
[13]. EDMs in particular provide quite stringent constraints on the low-energy CP- 
violating SUSY phases of the superparticle couplings. For certain values of the 
SUSY parameters, there is a disagreement with the experimental limits, resulting 
in the so-called SUSY CP problem. In this paper we make extensive use of the 
computer program CPsuperH2 . [TBI [TBI IT?] in choosing values for the various SUSY 
parameters in our processes. The most recent version of this program allows the 
user to compute EDMs, as well as various other low-energy observables. 

In a given decay, there are three types of CP-violating signals. First, there is 
the partial rate asymmetry. Any difference in the rate between process and CP- 
conjugate process is a signal of CP violation. As we will see, it is possible that 
the partial rate asymmetry is sizeable (of order 5-10% for the processes considered). 
Second, one has the modified (spin-dependent) rate asymmetry. Here one compares 
the rate for process and CP-conjugate process for the case in which the spin of 
one of the particles has been measured. This measurement is complicated, and so, 
given that the partial rate asymmetry can be significant, we do not consider the 
modified rate asymmetry in this paper. The third CP-violating signal is the triple- 
product (TP) asymmetry, which is proportional to v\ ■ {v% x £T 3 ) (each v% is a spin or 
momentum). This is due to terms of the form Tr[7 a 7 J g7 P 7 <J 75] in the square of the 
amplitude. Since only 3-body decays are considered here, a nonzero TP can arise 
only if a spin or polarization is measured. As noted above, such measurements are 
difficult, and so we do not consider TPs here. Thus, the only CP-violating signal 
analyzed in this paper is the partial rate asymmetry. 

All effects that violate CP require the interference of (at least) two amplitudes. 
The partial rate asymmetry is proportional to sin 5, where 5 is the relative strong 
(CP-even) phase between the interfering amplitudes. Strong phases can be generated 
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FIG. 1: Feynman diagrams for the decay x 2 — * Xi H\H\. In diagram (a) the decay is 
mediated by neutral Higgs bosons, while diagrams (b) and (c) are mediated by charginos. 



in one of two ways. First, one can have the exchange of gluons between the particles 
involved in the decay, leading to QCD-based strong phases. Unfortunately, we do 
not know how to calculate the strong phases in this case. Alternatively, the strong 
phases can be generated by the (known) widths of the intermediate particles in the 
decay. In the decays considered in this paper, the particles do not couple to gluons. 
Thus, the strong phase arises only due to the widths of the intermediate particles. 
This is good, given that we want to measure the CP violation, and not simply detect 
its presence. 

Now, SUSY includes two Higgs doublets which contain (in the gauge basis) two 
neutral scalars and one pseudoscalar. In the mass basis, these particles mix and one 
obtains three mass eigenstates Hi, H 2 and H 3 . In SUSY theories, the lightest mass 
is w,Hi = 0(100) GeV and we therefore take the final-state H to be Hi. 

Figure [l] shows the diagrams that contribute to yJi X\HiHi. The diagrams 
for the ZZ case are identical, but with Hi replaced by Z in the final state. The 
WW diagrams are similar, except that there is no diagram analogous to Fig. [T] (b) 
(we take "^4" to correspond to the W~ for this decay). Furthermore, the diagram 
analogous to Fig. [I] (c) involves intermediate neutralinos instead of charginos, and 
there is an extra diagram similar to Fig. [I] (a), but with the intermediate Higgs 
bosons replaced by the Z. Finally, for the ZHi case the diagrams are similar to 



Fig. 
Fig. 



1, but with one H\ replaced by Z. Also, there is an extra diagram similar to 
1 (a), but with the intermediate Higgs bosons replaced by the Z. Note that 



the neutralinos, x° (2 = 1, . . . 4), are the fermionic partners of the 7, Z, and neutral 
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Higgs bosons. 

The decay amplitudes that are of most interest to us arise from diagrams con- 
taining an intermediate particle that can go on-shell. Such diagrams can benefit 
from a resonant enhancement. The internal X121 H\ and Z can never be on-shell, 
while the H2 and H3 can be. (For a given set of SUSY parameters, the widths 
and T3 are calculable, as are the "off-diagonal widths" associated with transitions 
Hi Hj [18] . These terms are taken into account by CPsuperH2.0 when it com- 
putes the neutral Higgs propagator matrix.) The case of internal neutralinos (which 
appear in the case XY = W + W~) is more complicated. In principle, these could be 
on-shell. Unfortunately, we do not know their widths, making their contributions to 
the rate asymmetry uncertain. In practice, for the cases that we consider, the two 
heavier neutralinos have masses close to m~- and the two lighter neutralinos have 
masses close to or less than m~- . Thus, in the examples we consider, the neutralinos 
cannot go on shell and the uncertainty associated with the neutralinos' (unknown) 
widths is mitigated. 

In order for a given partial rate asymmetry to be appreciable, the interfering 
amplitudes should be of comparable sizes. Furthermore, since the rate asymmetries 
depend on an integration over phase space, the "large" contributions from these 
amplitudes should occur in similar regions of phase space. This latter requirement 
is met if the masses of the on-shell particles, H 2 and if 3, are similar. Fortunately, 
it is relatively common in SUSY theories that m# 2 ~ m# 3 [19j. The decay ampli- 
tudes of most interest to us, therefore, are those dominated by internal H2 and H3 
exchange. In calculating the square of the amplitude, however, we retain the con- 
tributions of various non- resonant diagrams (see Fig. [T]), since they can in principle 
give non-negligible contributions. With these ingredients, we calculate the partial 
rate asymmetry. Since it is assumed that the masses are known, we can deter- 
mine whether or not CP violation is likely in these decays. Measurement of these 
CP asymmetries will allow experimentalists to extract and/or constrain the SUSY 
parameters, including CP-violating SUSY phases. 

Analytical expressions for the amplitudes in question may be found in the Ap- 
pendix. Our notation is similar or identical to that employed in Ref. [15]. We have 
also used Feynman rules derived from Refs. [20l [21] where necessary. FeynCalc j'22.1 
was used to compute the squares of the amplitudes. The resulting expressions are 
quite messy and have not been included here. Please also note the following: 

1. Although only H2 and H3 can go on-shell, 1 we have allowed for the possibility 
of non-zero widths for other intermediate particles. For the purpose of our 
numerical work, all non-Higgs intermediate particles (with the exception of 
the Z, whose width is known) have been given a uniform width of 10 GeV. 



1 As noted above, the neutralinos can in principle go on-shell in the "WW" case, but they do not 
do so for the parameters that we use. 
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Since the intermediate particles in question are off-shell, the quantities that 
we compute should not be very sensitive to the value (s) assumed for the inter- 
mediate particles' widths. We have explicitly checked the effect of changing 
the uniform width from 10 GeV to 1 GeV and to 20 GeV for the data points 
plotted in Figs. [2] and [3] and have confirmed that the quantities shown in these 
plots are not very sensitive to such changes. 

2. The inclusion of intermediate particles' widths in our calculation implies 
that certain beyond-tree-level diagrams have effectively been taken into ac- 
count. Other such diagrams have not been included, leading to possible issues 
with gauge dependence or with respect to the CP-CPT connection [23| 124"] . 
CPsuperH2 . uses the Pinch Technique when computing the elements of the 
neutral Higgs propagator matrix [25l [26] . This technique was designed in such 
a way that certain classes of diagrams would be gauge independent. The neu- 
tral Higgs propagator matrix is a 4 x 4 matrix in this approach rather than 
the 3x3 matrix that one might expect. Rigorous application of the Pinch 
Technique for the present calculation is beyond the scope of our paper. In- 
stead, we have used the 3x3 physical Higgs block of the 4x4 propagator 
computed by CPsuperH2 . and have used Feynman rules consistent with the 
Unitary gauge. 

3. It is understood that Hi is not, in general, a CP-eigenstate, since the scalar 
and pseudoscalar Higgs bosons can mix when CP is violated. Nevertheless, we 
have checked numerically that the asymmetries for the H\Hi and ZH\ case 
are in fact zero when all CP-violating phases are zero. 

Having used the expressions in the Appendix to compute the widths for the 
processes and CP-conjugate processes, we subtract one from the other to obtain the 
respective partial rate asymmetries, 

a xy _ r fe~ -> XiXY) - T(xt -> xjXY) 
CP " r(x2 - XiXY) + T(xt - xtXY) ' [i) 

with XY = E X E X , ZZ, W + W~ or ZH X . 

The SUSY parameter space to be considered is enormous, and we have not 
made any attempt to perform a systematic parameter scan in our numerical work. 
Rather, we have chosen to focus on a small region of parameter space that is 
of interest. In particular, we have focused on a variation of the "CPX" sce- 
nario described in Ref. [17]. In the CPX scenario, the SUSY parameters /i, A t , 
A b and A T are chosen in such a way that \[iA t ^ T \ = 8M| USY , where Msusy 
represents the common mass scale of the third-generation squarks and sleptons 
(see Refs. [T7J I27J for further details), and A t;biT are trilinear couplings in the 
soft SUSY-breaking Lagrangian. This scenario was originally motivated by the 
observation that certain CP-violating terms in the neutral Higgs mass-squared 
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FIG. 2: Plots of partial rate asymmetries as a function of &A t - The dashed lines show 
the asymmetries for the case $i=$2=^ > 3=^A f) = ( J > A T = 0° (the values of the other SUSY 
parameters are given in the text). The scattered points indicate the asymmetries obtained 
by allowing all six phases to vary randomly between 0° and 360°. Low-energy experimental 
results have not been used to constrain the SUSY parameter space when generating these 
plots. 



matrix depend on lm(fiA t ) /M| USY . In our variation on the CPX scenario, we 
set \n\ = 0.6 TeV and Mq 3 =M c>3 =M^ 3 =M Z3 =M^ 3 =M SUS y=0.5 TeV, as well as 
|M 2 | = 2 | Mi | = 200 GeV, |M 3 | =\ TeV, m H ± = 300 GeV and tan/? = 5. Also, we 
set |A U) d )C , s ,e,/x|=|^4t,6,r| = 1 TeV and &A udcsefl = 0°. Our notation here is the same 
as that used in Ref. [TT] . 

Figure [2] shows the four asymmetries under consideration, plotted as functions 
of $A t - For the purpose of the plots, the hierarchy factors pj (J = Q,U, D, L, E) 
have all been set to 10. 2 The dashed line in each plot shows the result obtained by 



2 In the notation of CPsuperH2 . 0, the mass parameters Mj i and Mj 2 are assumed to be equal to 
each other (here J — Q,U, D, L, E and the subscripts "1" and "2" refer to the first and second 
generation, respectively). The third generation mass parameters are allowed to be different 
from those of the first two. The hierarchy factors are defined via the following expression: 

M .h, 2 =Pj M J 3 - 
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TABLE I: Constraints imposed when choosing SUSY parameter values for Figs. [3] and [4] 
The first five rows refer to the EDMs for Thallium, the electron, Mercury, the neutron 
and the muon, respectively; the sixth row contains the bound we enforce for the SUSY 
contribution to the muon anomalous magnetic dipole moment. Further discussion of some 
of the constraints may be found in the text. Note the following: (1) there is some variation 
in the confidence levels corresponding to the experimental upper bounds quoted in the 
second column; (2) Refs. [T7] and [28] contain further information regarding the Thallium 
EDM; and (3) to obtain the bound for the B — » X s j branching ratio we have combined 
errors in quadrature to obtain (3.52 ±0.25) x 10 -4 and have doubled the uncertainty. The 
range quoted in the table, and the constraint imposed, is thus at approximately the 2a 
level. 
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allowing $ J 4 t to vary between 0° and 360°, while keeping the other five phases set to 
zero. (I.e., we set $i=$ 2 =<3>3=<3>A i ,=3 ) A r = 0°, where $1,2,3 are the phases associated 
with the complex gaugino mass parameters Mi,2,3- Note that throughout this work 
we adopt the convention that $^ = 0°.) The scattered points in these plots show 
the asymmetries obtained by allowing all six phases to vary randomly between 0° 
and 360°. As is evident from the figure, for the parameters we have chosen, the 
asymmetries are strongly dependent on <3>^ t , although other phases contribute to 
the asymmetries as well. Analogous plots, showing the asymmetries as functions of 
the other five phases, do not demonstrate the same pronounced dependence on the 
other phases. 

In generating the data for Fig. [2] we have not made any attempt to impose the 
various low-energy experimental constraints that are available, since the purpose of 
the plots is to demonstrate the functional dependence of the asymmetries on the 
phases. We now turn to a more careful consideration of the asymmetries by also 
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taking into account several low-energy constraints. 

The constraints we impose are listed in Table |T| which also contains some com- 
ments regarding the constraints. We offer here a few additional comments. Let 
us first consider the muon anomalous magnetic moment. According to the au- 
thors of Ref. [33], the experimental value for the muon anomalous magnetic mo- 
ment exceeds the SM prediction by Aa M =a* xp - a™=(30.7 ± 8.2) x 1(T 10 , which 
represents a 3.7a deviation. One option would be to require that the SUSY con- 
tribution make up the difference between the experimental value and the SM pre- 
diction. We take a somewhat broader view and require the SUSY contribution 
to be between zero and (30.7 + 8.2) x 1(T 10 = 38.9 x 1(T 10 . The constraint 
quoted in Table [i is thus a^ USY = (19.45 ± 19.45) x 10" 10 . The last two rows 
of the table describe constraints on the SUSY contributions to AM^ and AMb s - 
The experimental values for these quantities are AMg^ = 0.507 ± 0.005 ps _1 and 
AMb s = 17.77 ± 0.12 ps _1 [35]. The corresponding constraints that we have listed 
are thus just the experimental uncertainties in these quantities. These constraints 
are tighter than they need to be, since we are ignoring the theoretical uncertainties 
within the SM. Nevertheless, these particular constraints are easily passed for the 
parameters we consider. One constraint that we have not directly imposed is on the 
ratio Rbtu = B(B~ — > t~v)/B sm (B~ — > t~v). There has been some discussion in 
the literature regarding the possible range of this ratio. For the parameters used to 
generate Fig. [3] (see below), we find Rb T u ~ 0.985, which is easily within the range 
derived, for example, in Ref. [36] . 

Figure [3] shows the four asymmetries plotted as functions of the mean width 
(defined to be [r(jfc -> Xi^Y) + T(xt -> XiXY)] /2, with XY = H X H X , ZZ, 
W + W~ or ZHi). The SUSY parameter values or ranges used to generate the plots 
in Fig. [3] are the same as those used to generate the scattered points in Fig. [2j 
but with three changes. First, instead of fixing the pj (J = Q, U, D, L, E) to a 
particular value, as was done for Fig. |2j we have now allowed the five hierarchy 
parameters to vary independently (and randomly) in the range 8 to 12. Second, we 
have imposed the constraints described in Table [Tj (The quantities in the table are 
computed automatically by the CPsuperH2.0 software.) Imposing the constraints 
from Table [I] leads to relatively strong limits on the SUSY parameter space. This 
led us to make a third change: to increase the efficiency of our numerical work, we 
allowed $2 only to take on the values 0° and 180°, and we restricted <E>A t so that it 
took on values between 140° and 220°. Of the 10 4 parameter sets that were originally 
generated in this manner, 432 were able to pass all of the cuts. The 432 parameter 
sets did not contain any cases in which $ 2 was 180° (i.e., only the parameter sets with 
$2 = 0° survived the cuts). Parameter sets that passed all of the constraints were 
used to compute the widths and asymmetries that appear in Fig. K3] 3 For the sets 



3 A few technical notes regarding our calculation are the following: (i) CPsuperH2 . follows three 
separate approaches when computing the neutron EDM and thus gives three separate estimates 
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FIG. 3: Scatter plots of partial rate asymmetries versus mean widths for X2 — ► XiXY, 
with XY = HtHx, ZZ, W + W~ or ZH X . 



of parameters considered, and respecting various low-energy bounds, asymmetries 
of order 10% are possible for \ 2 - > X\HiH\ and slightly smaller asymmetries occur 
for the other decay modes. 

With the parameters used for the plots in Fig. [3j the three neutral Higgs bosons 
had masses m Hl ~ 112.9-119.5 GeV, m H2 ~ 289.3-290.3 GeV and m H3 ~ 291.1- 
292.1 GeV. Also, the lighter and heavier charginos had masses m~ ~ 191.6 GeV 
and m~ 2 ~ 613.3 GeV, respectively. Thus, the parameters were such that the two 
heavier Higgs bosons could go on shell when mediating the various decay processes, 
as was noted above. 

Figure [4] shows the correlations between $ 3 and for the parameter sets that 
passed the constraints from Table [I] and that were subsequently used for the plots 
in Fig. [3j As is evident from the figure, there is a strong correlation between these 



for the EDM |17j . In our numerical work, we insisted that at least one of these three numbers 
satisfy the neutron EDM constraint listed in Table [I] (ii) We used the (approximate) default 
method for dealing with the dimension-six Weinberg operator, rather than supplying our own 
integration routine. See Ref. [17] for further discussion. 
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FIG. 4: Correlations between $3 and &A t - The points shown correspond to the same 
parameters used to generate Fig. |3j 

two phases that comes into play in allowing the constraints to be passed. 

It is useful to consider the observability of a 10% partial rate asymmetry. After 
considering the main contributions from the open two-body decay channels, we 
estimate the total width of the heavier chargino to be of order 10—20 GeV within the 
allowed parameter space region. This result, together with the \2 —* XiH\Hi partial 
width, allows us to obtain a representative value for the statistical significance of the 
CP asymmetry, S ~ \Acp\V2N \ where N is the number of events corresponding to 
the decay under study. (N represents the number of xt events, as well as the number 
of X2 events; these are assumed to be similar.) If, for example, the partial rate 
asymmetries are measured at an e + e~ linear collider, with a single \2 production 
cross section of order 20 fb and an integrated luminosity of C = 500 fb _1 jH [5], 
then the statistical significance turns out to be of order unity. This value should 
be understood as a very crude estimate, since the X2 production cross section is 
strongly dependent on the SUSY parameters and since we have not preformed a 
detailed analysis, nor considered the case of the LHC. Our point is simply to show 
that a 10% partial rate asymmetry could well be reachable. 

In conclusion, if SUSY is discovered in future experiments, it will become impor- 
tant to measure the underlying parameters of the theory. This paper has examined 
four chargino decay modes that could help provide insight into the CP nature of 
the theory. In particular, we have computed partial rate asymmetries for the decays 
xt -> X1H1H1, xt -> XiZZ, xf -> XiW + W~ and xt -> Xi ZH i. Rate asymme- 
tries have an advantage over some other CP-violating observables in that no spins 
or polarizations need to be measured. In the numerical example that we studied, it 
was found that the rate asymmetries for these decay modes are particularly sensitive 
to the phase of A t , although other phases contribute to the asymmetries as well. For 
the parameters considered, asymmetries of order 10% were found for xt ~^ X1H1H1; 
somewhat smaller asymmetries were found for the other decay modes. Should SUSY 
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be discovered, chargino decays could provide a useful avenue for investigating the 
CP structure of the theory. 
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APPENDIX: EXPRESSIONS FOR THE DECAY AMPLITUDES OF xt -» 

xfXY 

This Appendix contains analytical expressions for the amplitudes associated with 
the processes considered in this paper. 

Let us first clarify our notation for the various propagators. We define Breit- 
Wigner-type propagators for chargino and neutralino internal lines as follows, 

i(t + m)D(p\ m\ r) ee 2 ^ + ( A ' X ) 

pz _ m z _|_ 2I m 

where the tilde is used to distinguish the Breit-Wigner propagators from the 3x3 
Higgs propagator to be described below. We also employ the Breit-Wigner form 
of the propagator for the graphs mediated by a Z boson. The Z propagator in the 
unitary gauge is iD(p 2 z ,m 2 7,Tz)(—g al3 + PzPz I ' m D - 4 

As noted in the text, our calculation employs the 3x3 physical Higgs boson 
block of the full 4x4 neutral Higgs propagator computed by CPsuperH2 . 0. We also 
differ notationally from CPsuperH2 . in terms of the over-all normalization of the 
propagator. The specific relationship between the two sets of notation (within the 
physical 3x3 block) is the following, 

Dij(M 2 ) = Z)<f superH (M 2 )/M 2 . (A.2) 

Unless noted otherwise, our notation for coupling constants and diagonalization 
matrices follows the notation used in Ref. [15]. One exception is the definition of 



There has been some discussion in the literature regarding the correct form to use for the 
propagator of a spin-1 particle when the particle's width contributes to a rate asymmetry. See, 
for example, Ref. [37]. While this discussion is important in some contexts, we nevertheless use 
the "naive" form of the Z propagator in our numerical work, since the intermediate Z boson is 
far off-shell in the examples we consider. 
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9HiH z (which occurs in the u ZHi" decay), for which we use the notation defined in 
Ref. 0. 

In the following, we include explicit expressions for the \2 decays. The corre- 
sponding expressions for the xt decays can be obtained from the given expressions 
by taking the complex conjugates of the Lorentz-invariant pieces (B, C, D...), with 
the exception of keeping the propagator functions Dij and D unchanged. As an 
example, this procedure is demonstrated explicitly for the case of \2 —* XiHiHi- 



1. X,Y — H X H X 



The amplitude for x% "~ > Xi HxH% is given by 



M 



Hi Hi 



B + C7 5 ) + £{D + F 1 



MG + H^) u~ 2 (s llPl ) 



(A.3) 



where = (p 2 + p^Yi £ M = (p 2 + Pa)^ an d 

3 



B 



9 V"^ (-) 
z k=l 



D(£ 2 , m 2 ~ , T~ ) + Dip 2 , m 2 ~ , T- 



Xfc' Xfe' 



(A.4) 



C 



D 
F 
G 
H 



Z -J=T,9H^-D tJ (M 2 )g HjHlHl 



E^sW. D ^^L V xJ + D( P 2 ,m*~ ,r~ ) 



" fe=i 
z fe=l 



k=l 



fc=l 



(A.5) 



(A.6) 



fe=i 



-0 



/3 



\, :j nn,l \2 : !J 'i'h\: :j n ix ; X;j 



where we have defined M 2 = (p 3 +p 4 ) 2 , ry J k = g% ~+~-g c , 

p ' -"lXi Xfc t 

and the factors r/j are 771 = 6 and 772,3 = 2. Also, the m~ k denote the chargino 
masses. 
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The amplitude for xt ~^ XiH\Hi is given by 

M H ^ = u^(s 2 , P 2) [(B + C 7 5 ) + i (D + F7 5 ) + f> (G + Hr 



where 



(A.7) 



B 



C 



D 
F 
G 
H 



^ E .'/^ x : v . A, )•'/'/, /':'/ '/, 

V z i,_f=l 



„2 2 

9 \^_(-)« 



E%p' 



k=l 



D(£ 2 ,m 2 ~ ,r~) + D(p 2 ,m 2 ~ ,F 

v ^ Xfc ' Xfc y vr vi.' 



Xfc' Xfc' 



■7^ E ^x 2 -A,(M 2 )^^ lH] 

V ^ i,7=l 



fc=l 



Xfc' Xfc' 



Xfc' Xfc' 



z fe=i 
z fc=i 

-^E«B( P 2 , m i,r~), 

z k=l 

Z fc=l 



(A.8) 



Note that the spinors u~ i2 in Eqs. (A. 3) and (A.7) actually refer to u~- in both 
cases. Here, and throughout this work, we have manipulated expressions in such a 
way that the xt decay amplitudes are written in terms of the spinors u~- . 



2. A . Y = ZZ 



The amplitude for % 2 ~~ > Xi ZZ is given by 



M 



zz 



«xA S 2,P2j 



B + C^fT + 'ft'f (d + f 1 



+rh u {G + H 1 5 )+ri u {j + K 1 5 )}u~ 2 (s u p 1 )e x ^e x u >\ (A.9) 



where 



B 



g 2 m w 5 n n\A-2\ 

E 9 Hi z t z-D l3 {M )g H . 



vv 



(A.10) 



12 



g 2 



4 cos 2 e w £ 



=1 



c = 



ig 2 m w 



\/2 cos 2 9 W i j=1 ' 

,2 2 



~4cos 2 ^^ 



E™^^,^ 2 ,™- ,r~ ), 



D = 



"8 cos 2 e w ^ 

,2 2 



9 

= ^■o^0 u ^'-' ! - l - ,)iC ' 2 -" 



G = 



H = 



9 Z 



cos 2 9 W ^ 



E^, fc ^(p 2 ,^| fe ,r~ ), 



7 = y- ^^x^kfc 
^ 8cos 2 ^ 

^! 8cos 2 ^ 



^(^'<' r xJ-^(P 2 



2 -n 

' 1 Xfc^ 



and where we have defined oo^ k = UfyU^ ± UfyUfin, with 

^1) = + (cos 2 ^ - sin 2 9 w )S ij , 

% = iC R y\C R y^ + (cos 2 ^ - sin 2 6 W )^ . 



(A.H) 



^ 2 '<' r xJ-^ 2 '<' r xJ] > ( A - 12 ) 



(A.13) 



The unitary matrices Cl and C R are used to diagonalize the chargino mass matrix. 
The matrices U(l,r) appear in the chargino-chargino-Z Lagrangian as follows, 



9 



where P L{R) = (1 - (+) 7 5 )/2. 
3. X,Y = W+W~ 

The amplitude for x~2 — > x^W^W" is given by 



<"[£/*P L + t^PiilxJ^. (A- 14) 



WW 



~- u~(s 2 ,p 2 ) [(b + C 7 5 ) <T + 7^7" + F 7 5 ) + T'V (G + #7 5 ) 
+ (7^-7 M P3 + J (rfs -h)<r) (J + K 1 5 )]u~ 2 (s 1 ,p 1 )e^e x u ^ (A.15) 
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where P4 denotes the momentum of the W and where 



B = 

C = 

D = 

F = 

G = 



g 2 m w s n n\*i\ 

■ 2 3 

ig m w „ p 2 
— IT E 9 Bi xt% a^wD^M ) , 

VZ i,j=l 

2 E^,^(P 2 ^|o,r~o), 



fc=i 



,2 4 



5 



fc=l 

2 4 



E m x2^^(p 2 ? ^|o,r~o), 



2 ^ ' x k 
z fc=i 

z fc=i 



2 4 



J 



2 
2 



D(M 2 , m|, T z ) . 



(A.16) 



In the above expressions, m~o and T~o denote the neutralino masses and widths, 
respectively. Also, we define V^ k as V^ k = VfoVffi ± VfaVffi, where 



v/2 

The 4x4 unitary matrix N is used to diagonalize the neutralino mass matrix. 



V ( %) = (N) i2 (C R ) jl i=(N) u (C R ) j2 



4. X.Y = ZII\ 

The amplitude for X2 ~^ Xi ZHi is given by 

M ZH ' = u~(s 2 ,p 2 ) [(B + C7 5 ) p£ + (£> + F7 5 ) 7 M + (G + # 7 5 ) #7* 



(A.17) 



where p 4 denotes the momentum of the H 1 in the final state and where 



B = 



V2 cos 9w l= ^* H ^ 



E alz+z-DiAM 2 ) g Hl H jZ 
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2 o 

C = -p^ y 9u~+~-D ij (M 2 ) g H ,H-z 



"8co S 3^ w |^ Vy K+ m xJ 



^ 2 -^ 2 ]£(M 2 ,m|,r z ), 



2 g 2 m~ -D(p 2 , m~ , r v 

Xk ^ r ' Xfc 



£ = '"xfc^ '"S 

£r{ 4^2 cos W 

y ^ Xfc V s ' Xfc' Xk 

[c/i 2 + U% 



Uj k (g s M ~ + ~--igi ~ + ~-)+Ul k (g% ~+~-+ig* ~+~-) 

L K:JH lX k X 2 H ^Xl X 2 ' H Kt/H ixlx 2 ^ H ^xlx 2 ' 



g 2 m w 
+ ^ 5^ — Qh^vv 



8 cos 3 6*vi/ 
2 # 2 m~ Dip 2 , m 2 ~ ,T~ ) 

jfe 4^2 COS W 

v y Xfc vs ' Xk' Xk' 



D(M 2 ,m 2 z ,T z ), 



t) 



L yy HiXlx 2 y ^lXfeX 2 



x 2 y ^ix^x 2 y 



+ 



y ^ Xfc ^ ' > 

^ 4y^2 COS 9\y 

U 12 -U 12 ] D(M 2 ,m 2 z ,T z ), 



g 2 m w 



8 cos 3 6V 



2 , 2 D(e 2 ,m| fc) r~ ; 



^ = y^ 

jfe 4\/2cosfl w 

2 g 2 D(e,m 2 ~ ,T~ ' 
_ y ^ vs ' xfc' Xfc ' 

^ 4 v / 2cos^ 



U ^kxix k -<xtx-J + U n(9 S H ^- +i9H^\ > 
- <xix- k ) + U R ^kxix- k +<xtx~ k ) 



Xfc' X*' 



^ Ay/2 cos 9 w 
> 2 m 2 , r 



' Xfc' 



^ 4^2 cos 9 W 



U L(9kxix 2 +<xtx 2 ) + U n(9kxtx 2 -<xtx 2 \ ' 
U ^kxt^ + <xixi 



« Vy #iXfe X 2 y ^lXfc x 2 ; 



(A.18) 



As noted above, our notation for gn l H ] z is the same as that adopted in Ref. [2]. 
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